Study of low energy Sig and Cs implantation induced amorphisation effects in Si(lOO) 
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The damage growth and surface modifications in Si(fOO), induced by 25 keV Si^ cluster ions, 
as a function of fiuence, (p, has been studied using atomic force microscopy (AFM) and channehng 
Rutherford backscattering spectrometry (CRBS). CRBS resuhs indicate a nonhnear growth in dam- 
age from which it has been possible to get a threshold fiuence, 0o, for amorphisation as 2.5 x 10^^ 
ions cm~^. For cj> below (fio, a growth in damage as well as surface roughness has been observed. 
At a (/) of 1 X lO'^* ions cm~^, damage saturation coupled with a much reduced surface roughness 
has been found. In this case a power spectrum analysis of AFM data showed a significant drop, 
in spectral density, as compared to the same obtained for a fiuence, (j) < cj)o. This drop, together 
with damage saturation, can be correlated with a transition to a stress relaxed amorphous phase. 
Irradiation with similar mass Cs~ ions, at the same energy and fiuence, has been found to result in 
a reduced accumulation of defects in the near surface region leading to reduced surface features. 
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I. INTRODUCTION 



Cluster ion implantation can be regarded as a fore- 
runner technology as compared to the conventional ion 
implantation technique used to dope sub-micron devices 
Using cluster ions very shallow implantation can 
be achieved at very low energy. However, with cluster 
implantation, nonlinear effects arising in the energy loss 
processes, as a result of the correlated motion of the con- 
stituent atoms, play an important role in deciding the 
defect structure near the target surface. In addition to 
resulting in a nonlinear growth in subsurface damage, 
cluster ion impact, through sputtering, can also results 
in kinetic roughening and smoothening of the surface ex- 
posed [1]. In view of all this, there has been a lot of 
activities involving low energy cluster ion irradiation re- 
lated to nonlinear sputtering nonlinear damage and 
defect production [1, 0, 0] , along with the formation of 
various kind of surface features [8, 9, 10. fill, [l^. 

In connection with the above. Si, presents itself as a 
very important material where low energy cluster ions 
can be used for shallow implantation, of interest to tech- 
nology. In some earlier work, contrary to common ex- 
pectation, amorphisation upon ion irradiation has been 
shown to start from the surface rather than the ion pro- 
jected range [l^. Results of Molecular dynamics (MD) 
simulations with 5 keV Si, show that the ion impacts 
produce unrelaxed amorphous patches that have a fast 
quenched, liquid like structure [ij. With increase in 
ion fiuence these regions overlap producing a continuous 
amorphous layer [l5|. In fact, with increase in ion fiu- 
ence, there is a superlinear growth of amorphous volume 
fraction with a lot of stress build up in the matrix. At 



high fiuence there is an abrupt transition to a state with 
a flat amorphous-to-crystalline (a/c) interface [H, Il7| . 
In such a case, out of plane plastic flow with a reduction 
in the in-plane stress have been observed All this 

suggest that ion irradiation induced amorphisation in Si 
is more like a phase transition, initiated by a spontaneous 
collapse of the damaged region. Very recent MD simula- 
tions carried out by Marques et al show it to be initiated 
by a high concentration of interstitial-vacancy (IV) pairs 
or bond defects, formed in the system [l9|. Similar re- 
sults have also been shown by Nord et al zOj who have 
pointed out that the subsequent transition resulting in 
a uniform amorphous layer is neither a complete homo- 
geneous nor a complete heterogeneous mechanism. This 
makes Si an ideal system to study using low energy clus- 
ter ions where such a transition to a complete amorphous 
state is expected at a lower fiuence, primarily because of 
overlapping of collision cascades coming from constituent 
atoms. 

In the present paper we show some results of a sys- 
tematic study of the subsurface damage produced and 
the surface features generated in Si(lOO), from Si^ and a 
similar mass Cs~ ion implantation at 25 keV. Channel- 
ing Rutherford backscattering spectrometry (CRBS) and 
Atomic force microscopy (AFM) have been used for sam- 
ple characterization. Increase in cluster ion fiuence has 
been found to result in a nonlinear growth and saturation 
in damage leading to amorphisation. The transition to 
an amorphised state is found to be associated with a sig- 
nificant drop in the power spectral density of AFM data 
which initially increases with increase in fiuence. 



II. EXPERIMENT 
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Cleaned Si(lOO) wafers ( p-type, 1-2.5 Q.cm ) were ir- 
radiated with 25 keV singly charged negative ions viz 
Si^ and Cs" from a SNICS-II ion source (NEC, USA) 
using a low energy ion implanter facility. Mass analy- 
sis of the cluster ions was carried out using a 45° sector 
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FIG. 1: (Color online) (a) CRBS spectra for Si(lOO), implanted with 25 keV Sig" and Cs~ for different values of implantation 
fluence together with the same for a virgin sample, (b) Fluence {(j>) dependence of displaced Si atoms, Hd{<i>), for various 
implantations. Circles and triangles represent data for Sis and Cs implantations respectively. The continuous curve is a fit to 
the data points with clusters. 



magnet (ME/q'^ = 18 MeV amu). The base pressure in 
the target chamber during irradiations was maintained 
around 2x10"'' mbar. All the irradiations were carried 
out at room temperature with a beam flux of 2-3x10^° 
ions cm^^scc^^ (ion current of 2 — in A) at off the 
sample normal. In each case one part of the sample was 
kept unimplanted to serve as a reference. Five samples 
named S1-S5 were systematically irradiated with ions of 
similar mass (Si^ or Cs~) with gradually increasing ion 
fluence from 2x10^-^ cni^^ to 1x10^^ cm~^. Three of 
these, viz SI, S3 and S5 were irradiated using Si^ clus- 
ters to fluences of 2 x 10^^ cm~^, 1.2 x 10^^ cm~^ and 
1 X 10^^ cm~^ respectively. The remaining two samples, 
S2 and S4 were irradiated with 25 keV Cs^ ions to flu- 
ences of 1 X 10^^ cm~^ and 6 x 10^^ cm~'^ respectively. 
These data are shown in Table. 1. 

TABLE I: Sample names, ions used and integrated fluence. 

Samples fluence Ion species of 

25 keV total energy 



SI 


2 X 10'' cm" 


■i 


Sis" 


S2 


1 x 10'^ cm" 


2 


Cs- 


S3 


1.2 X 10'^ cm 


-2 


Sis" 


S4 


6 X 10" cm" 


2 


Cs" 


S5 


1 X 10" cm" 


2 


Sis" 



CRBS measurements were carried out on all the sam- 
ples with 1.35 MeV He+ with a Si surface barrier detector 
placed at 130° relative to the incident beam direction. 
The measurements were carried out at a steady beam 
current of 5 nA, using the 3 MV Pelletron accelerator 
(9SDH2, NEC, USA) facility at lOP, Bhubaneswar. In 
case of unirradiated Si(lOO), the reduction in the inte- 



grated total yield from random to a channeled spectrum 
was found to be ^5%. 

Following irradiation, the surface topography was ex- 
amined by AFM in the tapping mode, using a multi-mode 
scanning probe microscope (Nanoscope IV, Veeco, USA). 
Measurements were performed in ambient condition us- 
ing a Si cantilever with a nominal tip radius less than 
^ 10 nra. Image processing and analysis of the AFM 
data were carried out using the standard WSxM software 
package [21, 22]. 



III. RESULTS AND DISCUSSION 

The CRBS results as measured for all the flve samples 
viz. S1-S5 and a virgin sample (unirradiated area), are 
presented in Fig. 1(a). From the flgure, one can observe 
a gradual increase in surface peak intensity with increase 
in ion fluence, 0. This increase in the area of the surface 
peak, over and above that in the virgin sample, indicates 
growth in damage produced with increase in the number 
of displacements. In line with this, the lowest fluence ir- 
radiation is found to induce very little damage, too small 
to be seen through CRBS data. However, increase in ir- 
radiation fluence from SI to S3 (with cluster ions) and 
that for S2 to S4 (for Cs ions) are seen to result in en- 
hancements of the surface peak, indicating an increase in 
defect production near the surface. For the sample S5, 
with the highest cluster irradiation fluence, the surface 
peak is seen to be the most intense. This sample, with 
a total single atom fluence of 5 x lO""^"* cm~^ is expected 
to have an amorphous layer with a thickness of ^ 8nrn 
extending from the surface [2^. Cs^ irradiation to a 
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(a) 

♦ . 11.01 nm 




FIG. 2: (Color online) Surface morphology of samples (a) SI and 
of S2, I.e. 1 X 10^^ atoms cm"^ 



fluence, 0.6 times as that of S5, is seen to produced a 
damage distribution going deeper into the bulk. 

Now we look at the difference in the surface peak ar- 
eas between the irradiated and the unirradiated regions, 
which at lower irradiation fluence, is proportional to the 
number of displaced Si atoms, in the irradiated lattice. 
We denote this as Kd which has been estimated for vari- 
ous samples integrating data (fig 1 (a)) between channels 
340 and 440. Fig. 1(b) shows ion fluence, </>, dependence 
of Kd for various implantation fluence. From the figure, 
one can see a clear nonlinear growth and saturation in 
damage with increase in ion fluence. Similar results, re- 
garding a nonlinear growth in damage production and 
saturation has been seen earlier in case of 5 keV/atom 
Cn cluster irradiations of Si(lOO) [S]. 

To get an idea about damage saturation and the amor- 
phisation threshold, as has been done earlier 0, , we 
have fitted the three points as obtained for Sis irradiation 
to an equation of the form, 

= a{l ~ e~^^*°) (1) 

where a is a constant and (/)o corresponds to the threshold 
fluence for damage saturation. The best fit, as indicated 
by the smooth curve in Fig. 1(b), yields a value of 4>o 
equal to (2.5 ± 7%) x 10^^ cm~^. One can see, at a flu- 
ence of 6 X 10^'^ cm~^, the Kd value corresponding to 25 
keV Cs atoms is almost the same as that expected for a 
similar mass Sis cluster at the same energy. This means, 
the heavy single atom induced cascades produce almost 
the same amount of damage or defects as those generated 
by a similar mass cluster ion. At this high fluence, be- 
cause of overlapping of damage produced cluster induced 
nonlinear effects are difficult to be detected. However, in 
this case there are defects extending into the bulk while 



(b) 

2.72 nm 




(b) S2. Both samples have the same monomer fluence as that 



those for a cluster ion are better confined in a surface 
layer. 

The above (po value of 2.5 x 10^^ cm~^ as obtained 
for Si^ clusters, where damage saturation starts, corre- 
sponds to a total atomic fluence of ^ 1.25 x 10^^ cm~^. 
This agrees with the finding of Agarwal et al. [1^ who 
have shown the amorphisation threshold, for 5 keV Si in 
Si, to be 1 — 3 x 10^^ cm~^. In fact, with clusters, as 
expected the present value agrees with the lower limit. 
In view of this, at an Sis cluster fluence of 1 x 10^* cm~^ 
(corresponding to an atomic fluence of 5 x 10^^ cm~^), we 
are already well above the threshold for amorphisation. 



A. Surface features and AFM data 

As has been mentioned earlier, AFM has been used to 
study the surface topography in various samples. Some 
AFM pictures (top view), for (5 /xm x 5 /im) scanned 
areas, taken on the samples are shown in Figs. 2 and 
3. Figs 2(a) and (b) correspond to samples SI and S2, 
while Figs. 3(a) and (b) correspond to samples S3 and 
S4 respectively. 

Usually the features on irradiated surfaces are de- 
scribed through a height-height correlation function 
which contains three important roughness parameters: 
(i) the vertical correlation length cr, (ii) the lateral cor- 
relation len gth ^ and (iii) the roughness exponent a 
[H m, The lateral correlation length, ^, de- 

scribes the lateral characteristics of the surface, the 
roughness exponent a describing the static scaling prop- 
erties. The most commonly reported parameter of sur- 
face roughness i.e. cr, or the root-mean-square (rms) 
roughness, characterizes the surface, only along the ver- 
tical direction. This is defined as standard deviation of 
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FIG. 3: (Color online) Surface morphology of samples (a) S3 and (b) S4. Both samples have the same monomer fluence as that 
of S3, I.e. 6 X 10^^ atoms cm"^ 



the surface height profile, h(x,y), at each point {x,y) of times that of S3 is shown in Fig. 3(b). However, a much 
a reference surface plane from the mean height (< h >), smaller number of nanohillocks, compared to the S3 sam- 
as given by, pie, could be seen. It has a a value of 0.42nTO which is 

significantly smaller than that observed for S3. The en- 
hanced surface roughness as seen with Si^ implantations 
in SI as compared to S2 (and S3 as compared to S4), 
(2) even with lower fluence of irradiations, are primarily due 
to molecular effects coming from overlapping of collision 
cascades of constituent atoms. Such effects are absent 
where N is the^number of pixels, h, = h{x, y) being the during a single mass Cs irradiation which results in lower 
height at the « pixel. damage accumulation and lower surface roughness. Com- 

In case of the sample SI, as shown in Fig 2(a), one can pared to all the above, the a value of surface roughness 
clearly see black dots corresponding to nanometer sized observed for the sample S5, which had the highest clus- 
pits and bright spots corresponding to hillocks (color bars ter ion fluence, turned out to be 0.21nm. This was quite 
indicating the heights). The observed hillocks are seen gn^all compared to the same obtained for S3, irradiated 
to have heights ranging between 1 - 6 nm with average ^j^h an almost one order of magnitude lower fluence. 
height of 1.64nm. It has an rms roughness, a of 0.34 nm. However, the parameter a is not enough to give a full 

Compared to this, there was hardly any surface feature characterization of the surface because it is limited by its 
in the case of a pristine Si(lOO) sample (AFM data not sensitivity only in the vertical direction. For example, 
shown here). The rms roughness, g, of the pristine sam- images with exactly the same rms roughness values 

pie, was found to be 0.13 nm, which is about one third of ^an have different surface morphologies Hg. In view of 
that for SI, corresponding to a smooth polished surface. ^Yiis, a power spectral density (PSD) analysis is often 

used to look at surface features and their possible origin 
In case of the lowest fluence Cs implanted sample, S2, [s^, [Ml, HI, HI] . The PSD analysis is accomplished by ra- 
the a for surface roughness was found to be 0.21 nm, ly- dially averaging the square magnitude of the coefficients 
ing between that for SI (0.34nTO) and the pristine sample of the two-dimensional Fourier Transform of the digitized 
(0.13nm). What is important to see here is that Cs im- surface profile h{x, y) defined by 
plantation to a five times higher fluence is not able to 
generate a surface roughness as observed in case of SI. 

Compared to the above, S3, implanted with Sis clusters, 1 
to a fluence sixty times that of the sample SI, was found ~ 7^ 

to have an rms roughness, cr, of 1.21 nm. However, the 
above sixty- fold increase in Sis cluster fluence from SI where h{r) — h{x,y). Here q is the spatial frequency in 
to S3, did not result in a proportionate increase in the reciprocal space, is the scanned area of length L and 
number of nanohillocks in S3. The surface topography h{r) is the height at the position r. The PSD, which 
of the sample S4, implanted with Cs to a fluence five is the Fourier transform of the height-height correlation 
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exp < h{r) > 



(3) 
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FIG. 4: (Color online) (a) C(q), as a function of spatial frequency q, for Sis and Cs irradiated Si(lOO) at different irradiation 
fluence. Qi un-irradiated Si; V: SI (with Sis at 2 x lO^^cm"^); □, S2 (with Cs monomer fluence 1 x lO^^cm"^); A S3 (with 
Sis at 1.2 X 10"cm-2); 0, S4 (with Cs at 6 x 10^^ cm-^). (b) C{q) for samples S3 and S5. Inset shows the AFM image for S5 
for a scan size of 5/im x 5/im. 



function, thus turns out to be a function of only one 
parameter, spatial frequency, q. Identification of various 
processes for surface transport (s^l is carried out using 
the stochastic rate equations for the surface evolution 
[sst - In this analysis p^, C{q) = 2D{q) /^a^q'^], where 
D{q) is a term coming from noise correlation, a-y being 
expansion coefficients. The index 7 has values of 1,2,3 
and 4 representing four modes of surface transport viz. 
viscous flow, evaporation- condensation, volume diffusion 
and surface diffusion respectively. This index is further 
related to the roughness scaling exponent, a through the 
relation 7 = 2{a + 1). 

The PSD spectra, C{q), for all the samples SI to S4 
and pristine, are shown in Fig. 4(a). One can see, there 
is an increase in spectral strength in the order pristine, 
S2, SI, S4 and S3 where the pristine has the lowest and 
S3 has the highest values. This order is exactly the same 
if one orders the samples according to the a values for 
the surface roughness. S2 and S4 have five times higher 
ion fluence as compared to SI and S3 respectively. But 
the the surface features produced are seen to be much 
reduced. 

Having seen the above general features we now try to 
look at other details of C{q). One can see that the pris- 
tine sample {a = O.lSnm) has a cut off value, qc, which 
is about 0.001 nm^^. This corresponds to a large cor- 
relation length, ^ of 1 /iTO. Further, it has a 7 close to 
1. The Cs implanted sample, S2, has also a smooth sur- 
face (cr — 0.21nm). The corresponding C{q) indicates 
surface modulations over a similar length scale as the 
pristine sample. Compared to this, the Sis cluster im- 
planted SI sample, with one fifth of the fluence as in S2, 



shows modulations with a higher value of qc of the order 
of 0.006nm~^. This indicates a correlation length, of 
the order of 170nm. As shown earlier, it has a a value 
of 0.34nm. The Fourier index 7 is seen to be about 2.5. 
This yields an a value of 0.25, indicating the surface to 
be self affine with anisotropic scaling along lateral and 
perpendicular directions. 

The C{q) for the higher fluence Cs irradiated sample, 
S4, shows a qc which is almost similar to that for the lower 
fluence Sis implanted sample, SI, indicating a similar 
correlation length, ^, of 170nm. However, it shows a 7 
value of around 2.2 resulting in an a value ~ 0.1. With a 
a value of 0A2nm, it has a higher roughness. This means, 
with a 300 times higher fluence, Cs ions generate similar 
surface features as obtained for a low dose Sis implanted 
sample, SI. Compared to this, the Sis cluster implanted 
sample S3, with a fluence which is one fifth of that in 
S4, shows a 7 value close to 4, indicating an a value 
close to unity. This indicates the surface modulations to 
be self similar. But this sample has the highest surface 
roughness, a of 1.21nm with an average height of 5nm. 

Now we look at what happens when cluster fluence 
is increased to a value well beyond the amorphisation 
threshold as in case of S5. An AFM image of the top 
view of S5, taken with a (5fim x 5fim) scan size is shown 
in Fig. 4(b). The C{q) spectrum for S5 is also shown 
along with that for S3 in the same figure. One can clearly 
see that increasing the cluster ion fluence from 1.2x10^^ 
to IxlO'^* cm~^, in going from S3 to S5, has resulted in 
no further change in the 7 value which has been found 
to saturate at 4. But the surface modulations changed 
to have a a value of 0.21nm with a correlation length 
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^ 125nm. With an a value close to 1, the surface 
features in S5 are self similar indicating isotropic scal- 
ing. However, as compared to the S3 case, C{q) shows a 
significant reduction in the magnitude. At qc, the ratio 
between the two is about 33. This can also be seen using 
the formula C{q) — [aa^£_^/n] at Qc — !/£,. Since S3 and 
S5 have almost same a as well as ^ values, the ratio of 
the C(g)s turn out to be nearly the same as the ratio of 
the for the two cases. This is seen to be (1.21/0.21)^ 
which is just about right. It is important to mention 
that the surface modulations in S5 show a mean height 
of O.Gnm. The small value of C{q) for S5, is therefore 
seen to be coming from a correlation between small but 
nearly equal heights with a small a value. Compared to 
this C{q) in S3 comes from a correlation between higher 
heights, with an average value of 5nm, again with a much 
higher value of a. 

Earlier, 5 keV Si impact on Si has been shown to re- 
sult in creation of amorphous poc kets coming form local 
melting and rapid quenching The stress produced 

can result in formation of pits and bumps on the surface. 
This is in addition to roughening resulting from sputter- 
ing. This local melting and the associated movement of 
atoms, at lower fluence as in SI, may be responsible for a 
7 value between 2 and 3. It also leads to a higher surface 
roughness as compared to a pristine sample. At lower flu- 
ence there are sparsely distributed amorphous pockets in 
the matrix. Increase in the implantation flucncc results in 
a fast growth in the number of amorphous patches result- 
ing in a growth in surface roughness, a. This also results 
in a growth in mean height, < h >, of surface structures 
produced, resulting in a growth of height-height correla- 
tion. Merging of amorphised regions at higher fluence, 
results in a building up of stress from a large number of 
bond defects. Finally there could be stress relaxation as 
the damaged lattice becomes unstable. This could result 
in a transition to a state with smaller surface features 
which is achieved by an effective movement of atoms in a 
lateral direction. This could be the reason behind getting 
a 7 value close to 4 as proposed for surface diffusion. This 
is probably how a smooth amorphous to crystalline in- 
terface can occur in Si under high fluence ion irradiation 
[TtI . This way crystalline to amorphous transition in 
Si, upon ion irradiation, is more like a phase transition 
induced by an accumulation of sufficient number of de- 
fects which was also suggested by several groups earlier 



0, In the present case, onset of this occurs at 

a cluster fluence of around 2.5x10^'^ cm~^. For a much 
higher cluster fluence (as in S5) , a continuous amorphous 
layer parallel to the surface is produced, leading to much 
reduced values for a and C{q). It is therefore not sur- 
prising that the ^, a and 7 values as obtained for S5 are 
very similar to those obtained for a-Si films [2^. This 
also confirms that the top surface of the S5 sample is 
actually amorphised, in agreement with channeling data 
(Fig. 1(b)) where saturation in damage production has 
been obtained. It is also important to realize that for 
S4 is much higher than that in case of S3 (Fig. 1(b)). In 
fact it is closer to that in S5. However, complete amorphi- 
sation does not occur here because the defects produced 
by Cs implantation are distributed over a greater depth 
resulting in comparatively less defect accumulation near 
the surface. 



IV. CONCLUSION 

To conclude, we have carried out a systematic study 
of 25 keV Si^ implantation induced damage and surface 
modifications in Si(lOO) where a nonlinear growth in sub- 
surface damage, with fluence, is observed. The damage 
produced by similar mass Cs^ ions, of the same energy, 
is seen to be distributed over a greater depth leading to 
much reduced surface features. With Sis clusters, the 
threshold fluence for amorphisation of Si surface is found 
to be 2.5 X 10^^ clusters cm~^, in agreement with ear- 
lier published data. Most importantly, at higher cluster 
fluence a transition to an amorphous state resulting in a 
much reduced surface roughness is indicated. 



V. ACKNOWLEDGMENTS 

The authors would like to thank A.K. Behera for the 
efficient running of the implanter facility and all the op- 
erators of accelerator laboratory, lOP for help during 
CRBS runs. The technical help regarding some prelim- 
inary AFM measurements and checks by S. Varma of 
lOP and S. Bhattacharjee of IMMT, Bhubaneswar are 
also gratefully acknowledged. 



[1] V. N. Popok and E. E. B. Campbell, Rev. Adv. Mater. 

Sci. 11, 19 (2006). 
[2] I. Yamada, J. Matsuo, N. Toyoda, and A. Kirkpatrick, 

Mater. Sci. Eng. R: Reports, 34, 231 (2001). 
[3] Y. Teranishi, K. Kondou, Y. Fujiwara, H. Nonaka, M. 

Tomita, K. Yamamoto, T. Fujimoto, and S. Ichimura, 

Jpn. J. Appl. Phys. 45, 5528 (2006). 
[4] H. H. Andersen, A. Brunelle, S. Della-Negra, J. Depauw, 

D. Jacquet, Y. Le Beyec, J. Chaumont, and H. Bernas, 



Phys. Rev. Lett. 80, 5433 (1998); A. V. Samartsev, A. 
Duvenbeck, and A. Wucher, Phys. Rev. B 72, 115417 
(2005). 

[5] L. Shao, M. Nastasi, X. Wang, J. Liu, and Wei-Kan Chu, 

Nucl. Instrum. Methods Phys. Res. B 242, 503 (2006); 

J. Liu, X. Wang, L. Shao, X. Lu, and W. K. Chu Nucl. 

Instrum. Methods Phys. Res. B 190, 787 (2002). 
[6] H. P. Lenka, B. Joseph, P. K. Kuiri, G. Sahu, and D. P. 

Mahapatra, Nucl. Instrum. Methods Phys. Res. B 256, 



7 



665 (2007). 

A. I. Titov, A. Yu. Azarov, L. M. Nikulina, and S. O. [23 
Kucheyov Phys. Rev. B 73, 064111 (2006). 
S. Prasalovich, V. N. Popok, P. Persson, and E. E. B. 
Campbell, Eur. Phys. J. D 36, 79 (2005). [24 
J. Samela, K. Nordlund, V. N. Popok, and E. E. B. 
Campbell, Phys. Rev. B 77, 075309 (2008). 
H. S. Park, H. J. Jung, and W. K. Choi, Thin Solid Films [25 
475 36 (2005). 

L. P. Allen, D. B. Fenner, C. Santeufemio, W. Brooks, [26 
and I. Yamada, J. Appl. Phys. 92, 3671 (2002). 
V. N. Popok, S. V. Prasalovich, and E. E. B. Campbell, [27 
Nucl. Instrum. Methods Phys. Res. B 207, 145 (2003). [28^ 
K. Nakajima, H. Toyofuku, and K. Kimura, Jpn. J. Appl. 
Phys. 40, 2119 (2001). [29 
T. Diaz dela Rubia and G.H. Gilmer, Phys. Rev. Lett. 
74(13), 2507 (1995). [30 
S.U. Campisano, S. Coffa, V. Raineri, F. Priolo, and E. 
Rimini, Nucl. Instrum. Methods Phys. Res. B 80/81, [31 
514 (1993). 

O.W. Holland, S.J. Pennycook, and G.L. Albert, Appl. [32 
Phys. Lett. 55, 2503 (1989). 

G. Bai and M.A. Nicolet, J. Appl. Phys. 70, 649 (1991). [33 
C. A. Volkert, J. Appl. Phys. 70, 3521 (1991). 
L. A. Marques, L. Pelaz, P. Lopez, I. Santos, and M. [34 
Aboy, Phys. Rev. B 76 153201 (2007); L.A. Marques, [35 
L. Pelaz, M. Aboy, L. Enn'quez, and J. BarboUa, Phys. 
Rev. Lett. 91, 135504 (2003); L.A.Marques, L. Pelaz, J. [36 
Hernandez, J. BarboUa, and G.H. Gilmer Phys. Rev. B 
64, 045214 (2001). [37 
[20] J. Nord, K. Nordlund and J. Keinonen, Phys. Rev. B 65, 

165329 (2002). [38; 
[21] WSxM free software. Avilable at www.nanotech.es 
[22] 1. Horcas, R. Fernandez, J.M. Gomcz-Rodrigucz, J. 
Colchero, J. Gomez-Herrero, and A. M. Baro, Rev. Sci. 



Instrum. 78, 013705 (2007). 

A. Agarwal, T. E. Haynes, D. J. Eaglesham, H.J. Goss- 
mann, D.C. Jacobson, J. M. Poate, and Y. E. Erokhin, 
Appl. Phys. Lett. 70, 3332 (1997). 

M. DobeU, R,M. Ender, U.S. Fischer, M. Suter, H.A. 
Synal, and D. Vetterli Nucl. Instrum. Methods Phys. Res. 
B 94, 388 (1994). 

E.A. Eklund, E.J. Synder, and R.S. Williums, 6285, 157 

(1993). 

K.R. Bray and G.N. Parsons, Phys. Rev. B 65 035311 
(2001). 

D.B. Fenner, J. Appl. Phys. 95, 5408 (2004). 

S.J. Fang, S. Haplepete, W. Chen, C.R. Helms, and H. 

Edwards, J. Appl. Phys. 82, 5891 (1997). 

S.J. Fang, W. Chen, T. Yamanaka, and C.R. Heims, 

Appl. Phys. Lett. 68, 2837 (1996). 

R. Petri, P. Brault, O. Vatel, D. Henry, E. Andre, P. 
Dumas, and F. Salvan, J. Appl. Phys. 75, 7498 (1994). 
W.M. Tong, R.S. Williams, A. Yanase, Y. Segawa, and 
M.S. Anderson , Phys. Rev. Lett. 72, 3374 (1992). 
C. Ruppc and A. Duparre; Thin Solid Films 288, 8 
(1996). 

M. Senthilkumar, N.K. Sahoo, S. Thakur, and R.B. 

Tokas, Appl. Surf. Sci. 252, 1608 (2005). 

C. Herring, J. Appl. Phys. 70, 3521 (1950). 

W.M. Tong, E.J. Synder, R.S. Williams, A. Yanase, Y. 

Segawa, and M.S. Anderson, fi277, L63 (1992). 

S.G. Mayr and R.S. Averback Phys. Rev. Lett. 87, 

196106 (2001). 

M.L. Swanson, J. R. Parsons and C.W. Hoelke, Radiat. 
Eff. 9,249 (1971). 

T. Motooka, S. Harada, and M. Ishimaru, Phys. Rev. 
Lett. 78, 2980 (1997). 



